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A fully digital data acquisition system based on a field-programmable gate array (FPGA) was developed for 
a CsI(TI) array at the External Target Facility (ETF)in the Heavy Ion Research Facility in Lanzhou(HIRFL). To 
process the CsI(TI) signals generated by y-rays and light-charged ions, a scheme for digital pulse processing 
algorithms is proposed. Every step in the algorithms was benchmarked using standard ~y and œ sources. The 
scheme, which included a moving average filter, baseline restoration, leading-edge discrimination, moving win- 
dow deconvolution and digital charge comparison was subsequently implemented on the FPGA. A good energy 
resolution of 5.7% for 1.33 MeV y rays and excellent a-y identification using the digital charge comparison 
method were achieved, which satisfies CsI(Tl) array performance requirements. 
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I. INTRODUCTION 


The structure of atomic nuclei near drip lines is one of the 
most fascinating fields for nuclear physicists, and has con- 
tinuously attracted researchers to build large facilities for ex- 
perimental studies [1, 2]. One such facility is the External 
Target Facility(ETF) at the Heavy Ion Research Facility in 
Lanzhou(HIRFL) [3]. The ETF is a large integrated exper- 
imental platform for nuclear physics research. It comprises 
several detector systems and can provide complete kinematic 
measurements of nuclear reactions at intermediate and high 
energies [3-7]. The CsI(T]) array, consisting of 1024 CsI(T) 
detectors, is one of the most important detector systems in 
the ETF. It can measure y-rays up to 10 MeV in the center- 
of-mass(CM) frame [8]. Such a design, an inorganic scintil- 
lation array with high granularity, is also used in many detec- 
tor arrays, such as DALI2 at RIKEN [9, 10] and CALIFA at 
FAIR [11, 12], and it enables a relatively good energy reso- 
lution of the y-rays to be obtained owing to its good angular 
resolution [10]. 

In combination with silicon strip detectors, the CsI(TI) ar- 
ray at the ETF can also be used to measure light-charged 
particles. However, the energy loss of light-charged parti- 
cles in CsI(Tl) crystals is several tens of times greater than 
that of y-rays. To satisfy the highly dynamic requirements of 
the measurements, the electronic system was updated. High- 
granularity detectors require highly integrated electronic sys- 
tems [13-16]. The traditional solution, also known as 
the previous scheme of the CsI(TI) array, uses application- 
specific integrated circuit (ASIC) chips for signal process- 
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ing [17]. However, these chips are highly customized and it is 
difficult to extend their functionality. Waveform digitization 
based on flash analog-to-digital converters (ADC) has been 
developed for over two decades. The analog signals extracted 
from the detector or charge sensitive amplifier (CSA) are di- 
rectly digitized. Therefore, a minimal loss of signal infor- 
mation can be achieved using a much simpler circuit. Other 
advantages include a higher sustained count rate, flexibility 
through a variety of digital pulse process (DPP) algorithms, 
and compact structures, which facilitate the development of a 
highly integrated system. Owing to these aforementioned ad- 
vantages, a fully digital technique was chosen as the solution 
for a new measurement system for the CsI(TI) array. 


However, this scheme has one drawback. The digitization 
of the entire waveform means that a significant amount of 
data must be read. For a one-channel signal digitized by a 
14-bit flash ADC with a sampling frequency of 50 MS/s, the 
data rate was approximately 83.4 MByte/s. When this flash 
ADC is used for the CsI(T]) array at the ETF, the amount of 
data is extremely large, and the transmission bandwidth of 
the Data Acquisition (DAQ) system may be overloaded. Zero 
suppression should be performed after wave digitization to 
reduce data transmission pressure. Using DPP algorithms in 
the onboard field-programmable gate array (FPGA) to extract 
specific waveform information, such as the amplitude and ar- 
rival time, the entire amount of waveform sampling data can 
be reduced to a few physical quantities, further reducing the 
data volume. Moreover, owing to conflicts between the FPGA 
limited computational resources and the large number of de- 
tector channels, the new electronic system may require a com- 
promise in precision with somewhat uncomplicated DPP al- 
gorithms, and these algorithms are the focus of this study. 

The remainder of this paper is organized as follows: 
Section II concentrates on general considerations regarding 
which procedures should be executed in the FPGA and in 
what order. Section II lists the specific DPP algorithms for 
each procedure and discusses the rationale for algorithm se- 


67 lection through bench testing. Section IV presents the final 
e DPP algorithms scheme in the FPGA and discusses their per- 
e formance after implementing all the algorithms in the FPGA. 
zx Section V summarizes the study findings with an outline of 
71 future research directions. 


7z Il. GENERAL CONSIDERATION OF DPP ALGORITHMS 


7s In nuclear physics, signals from detectors are completely 
74 random. This indicates that DPPs for the CsI(T]) array should 
75 perform well in the time domain [18]. As a calorimeter, the 
76 critical physical quantity measured by the CsI(TI) array is the 
77 energy of the incident gamma rays or light ions. Thus, the 
7s method for obtaining energy information from the recorded 
79 waveform of the detector is the most critical aspect of this 
so Study, and energy resolution is a key criterion for DPP algo- 
s1 rithms. Further algorithms such as a smooth filter and base- 
s2 line recovery were used to achieve a good Signal to Noise 
83 Ratio (SNR) value and energy resolution. The signal arrival- 
s4 time information is responsible for generating a system trig- 
es ger. It also helps to reduce the background noise by selecting 
sæ appropriate time windows for off-line data analysis [19]. 
s7  CsI(Tl) crystal has been used for many years to identify 
, es light ions and gamma rays using pulse-shape analysis. This is 
s9 because the responses of the CsI(T1) crystal vary with differ- 
æ ent types of particles, resulting in different waveform shapes 
ə being generated by the CsI(T1) detector [20-24]. In the case 
ə2 Of a CsI(TI) array with high spatial coverage and large gran- 
ə ularity, the gamma rays and light ions can strike different el- 
ə ements of the detector simultaneously and form different hit 
ə Clusters. When the energy spectrum is reconstructed by some 
æ algorithms, such as the add-back technique, the spectrum be- 
97 tween the charged particles and the gamma rays may be su- 
ı 9s perimposed. Thus, clear separation of the gamma rays and 
» 99 charged particles can result in a good gamma energy spectrum 
with a lower background of charged particles. Pulse-shape 
analysis can improve the performance of particle identifica- 
tion(PID) compared to the traditional AE-E method [25]. 
In some situations, when the atomic number Z of the mea- 
sured ions is less than four, pulse shape analysis can also 
help simplify the detector setup because the CsI(T1) detec- 
tor can perform identification by itself [26-28]. Considering 
the aforementioned advantages, on-line algorithms for pulse 


shape analysis are also in demand. 
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Fig. 1. DPPs in FPGA for the CsI(TI) array of ETF. 
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109 Overall, In this study, DPP algorithms are organized in the 
110 FPGA in the manner shown in Figure 1. The final outputs 


111 were the arrival time, energy, and quantities representing the 
incident particles PID results. Each process as shown in Fig- 
ure | is described in detail in the following sections. 
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HI. ALGORITHMS SELECTIONS WITH OFF-LINE 


ANALYSIS 


To identify suitable algorithms for each process shown in 
Figure 1, a test bench including an element of CsI(T1) array 
and the DAQ system was set up. The DAQ system is an inter- 
mediate development product of the CsI(T1) array that con- 
tains a CSA module and a DAQ board. A 14 bit flash ADC 
with a sampling rate of 50 MS/s was embedded in the DAQ 
board to digitize the signal from the CSA module. The DAQ 
board can be operated in two modes: raw waveform mode, in 
which the Flash ADC data are recorded directly into PC mem- 
ory, and algorithm mode, in which the data are processed by 
algorithms in FPGA with only the results recorded. Because 
the DPP algorithms are not specified, the test bench is oper- 
ated in the raw waveform mode to select the appropriate algo- 
rithms for the CsI(T]) array. The criteria for selecting suitable 
algorithms are based on the tradeoff between performance, 
FPGA source consumption, and execution speed. The Flash 
ADC sampling rate was set to 25 Ms/s in the FPGA because 
the data volume was too large for the raw waveform mode. 


A. Smooth filter 


Because the detector signals are always distorted by ran- 

dom noise, the goal of this procedure is to reduce the high- 
frequency noise while not significantly changing the detector 
signals to preserve the signal characteristics and improve the 
data SNR. This procedure was performed initially because it 
may slightly alter the raw waveform and affect the perfor- 
mance of the other procedures. 
142 The moving average filter (MAF) is the most commonly 
used smooth filter in the time domain owing to its sim- 
144 plicity, ease of implementation, and rapid execution speed. 
145 In addition to the aforementioned advantages, it offers the 
lowest noise for a given edge sharpness for any linear fil- 
ter [18]. Other optional filters include the Savitzky—Golay 
[29], binomial [28, 30], Whittaker [31], and the Kalman fil- 
ters [32, 33]. Although these smooth filters perform well in 
multiple fields such as Raman and Mossbauer spectroscopy, 
complex algorithms make them difficult to implement in the 
FPGA. For this reason, the MAF was chosen as a smooth fil- 
ter for the CsI(TI) array. 

The only parameter for MAF is the number of samples av- 
eraged over. Its value should be chosen with care because as 
the value increases, the noise decreases, while the edges be- 
come less sharp. The aforementioned test bench, which was 
operated using cosmic rays, was used to acquire raw wave- 
159 form data, and the MAF algorithms were executed with dif- 
160 ferent parameter values in the off-line analysis. The results 
161 are presented in Figure 2. The reason for using a power of 
162 two for the parameter values is that the division is simpler to 


163 implement in the FPGA using only shift operations. Finally, 19 rectly. Therefore, both methods are treated as a single method 
164 Owing to its good performance in high-frequency noise reduc- 200 and are hereafter called iterative methods. 
165 tion and small waveform changes, the value of the parameter 
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Fig. 2. (Color online) Comparison among MAF results with differ- 
ent values of N. Inset illustrates leading edges in the picture. 


B. Baseline restoration 


In many cases, the baseline is assumed constant at all times. 
This enabled the baseline value to be measured at any time 
while the DAQ system was idle during the experiment. How- 
ever, this assumption is only true when the signal length is 
short and the count rate is relatively low. This is certainly 
not the case in this instance, because the CSA signal ex- 
tracted has a long tail. A good solution is to use optimal 
filters [34]. However, complicated arithmetic discourages us 
from proceeding further. To simplify the procedure, it was 
assumed that the baseline remained constant within each se- 
lected dataset, which contained only one signal waveform 
with several baseline data points in front. This baseline data 
can subsequently be used to calculate the baseline value of 
each individual event. Reference [35] showed that two values, 
the average and median of the baseline data, can be used to 
evaluate the baseline level. The method for obtaining average 
baseline data is clear. The median is the exact quantity in the 
middle of the baseline dataset when ordered. Reference [35] 
showed that the median is a better estimate of the baseline 
level than the average over a wide range of count rate loads. 
The other two methods, which are called "averaging over the 
selection set" and "averaging over the flat chunk selection set" 
respectively, are also introduced in [35]. These two methods 
are identical in terms of how to proceed, albeit differ in terms 
of data selection. Further detail regarding both methods can 
be found in [35], and the main procedures are summarized 
in the flowchart shown in Figure 3. Over a wide range of 
count rate loads, these two methods provide better estimates 
of the baseline than either the average or the median alone. 
In this case, because the baseline dataset has already been 
selected, the procedure list in Figure 3 can be processed di- 
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Fig. 3. The main procedures of the methods, averaging over the 
selection set and averaging over the flat chunk selection set, from 
[35]. 


Here, there are three methods: the average of the base- 
line dataset, median of the baseline dataset, and the iterative 
method. To evaluate performance, the baseline samples in the 
data used in the previous section were processed using these 
methods. The results of the baseline recoveries are shown 
in Figure 4. The discrete and continuous histograms are dis- 
played using the median and average methods, respectively, 
because the data types for these two methods are "int" and 
“float.” The histogram obtained using the iterative method 
is similar to that obtained using the median method. The 
only difference is the width of each discrete part, which is 
equal to twice the minimum difference parameter set in the 
algorithm(é, in Figure 3, here the value is 0.1, further details 
can be found in [35]). Gaussian functions were used to fit the 
envelopes of the three histograms; the results are listed in Ta- 
ble 1. Because the data are all integers for the median method, 
the most probable value(MPV) of the histogram is treated as 
zero and not the mean fitting parameter. From a comparison 
of the MPVs of the three methods, it can be concluded that all 
three methods are good estimators of the baseline levels, and 
therefore, fully meet the performance requirements 

Algorithms subsequently become the focus of the selection 
criteria. The iterative method was the most complex of the 
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three methods. However, this method is powerful because 
it can continue to provide reasonable baseline values when 
the signal samples are also included the dataset [35]. The 
median method is simpler, albeit the data ranking algorithm 
is not "FPGA friendly.” Therefore, for the CsI(Tl) array DAQ 
system, the average method is preferred. 
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Fig. 4. Typical results of baseline restorations obtained with the fol- 
lowing algorithms: (a) Median, (b) Average, and (c) Iteration. 
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TABLE 1. Fitting parameters of the three methods. 


Methods MPV Sigma 
Median 0 3.175 
Average 0.0008 3.147 
Iteration 0.007 3.187 


C. Arrival time 


According to the ETF design scheme, the CsI(TI) array was 
intended to generate a trigger for the entire system. With the 
complete digitization of the input waveform and the powerful 
calculation capability of the FPGA, the CsI(TI) array trigger 
signal can be generated and controlled by software. This sim- 
plifies the electronic system by eliminating the necessity for 
additional electronics, such as splitters and time discrimina- 
tors. 


An ETF trigger system consists of 2-level trigger genera- 
tors [36]. At the front end, primary trigger signals were gen- 
erated according to the signal shapes and logical relationships 
between the readout elements of a particular detector included 
in the trigger system. These primary trigger signals were sub- 
sequently fed into the global trigger logical unit to generate 
an event trigger signal for the entire system according to the 
physical interests. To improve logic operation effectiveness 
for signals with different delay times and time jitter, each pri- 
mary trigger signal can be delayed and widened separately 
in a global trigger logical unit. From this perspective, time 
resolution is not the key element for the CsI(TI) array. 


Several methods are available to determine the arrival time 
of a signal. Common methods, which are identical to those 
in the analog scheme, are leading-edge discrimination (LED) 
and constant fraction discrimination(CFD). LED is the sim- 
plest method; however, it has a large time jitter owing to 
the time-walk effect. To achieve better performance, cor- 
rection should be performed, which is a significant task for 
the FPGA programmer. Therefore, CFD was introduced to 
eliminate the time walk effect. There are two methods for 
implementing this algorithm in an FPGA: constant-fraction 
zero-crossing (CFDzc) and digital contact-fraction discrimi- 
nation (dCFD) [37]. CFDzc is a digital version of the classic 
analog CFD [38], whereas dCFD is similar to an LED with a 
different threshold value, equal to the constant fraction of the 
signal amplitude [39]. Better performance in terms of time 
resolution was obtained with the CFDzc method than with 
the dCFD method [40]. Other methods exist for determining 
the arrival time of a signal, such as the RC-CR2 filter [41] and 
pulse-shape fitting. However, utilization of these algorithms 
is resource-intensive in the FPGA and considered outside the 
scope of this research. 

Because time resolution is not the key criterion for the 
CsI(TI) array, the LED method without interpolation was cho- 
sen to determine the arrival time, considering the expected 
high computational resource consumption of the energy ex- 
traction and PID procedures. Figure 5 shows the relationship 
between the amplitude of the waveform and arrival time using 
the previous data. The arrival time is the difference between 
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the reference time and the time when the first sample point 
is above the threshold. The time reference was selected as 
the point at which the ADC waveform passed through 90% 
of its full amplitude at the leading edge. Linear interpolation 
was used to reduce the reference time jitter. A clear depen- 
dency between the waveform amplitude and the arrival time 
is shown in Figure 5. Another phenomenon is that the lower 
the waveform amplitude , the greater the time jitter. This is 
because a low amplitude causes the leading edge of the wave- 
form to become flat, and it is difficult to determine the time at 
which the waveform crosses the threshold and the time refer- 
ence. Although the time jitter is as large as approximately 
3 us when the waveform amplitudes are extremely small, 
this may not be useful because a threshold can be set in the 
FPGA to determine the recorded signals. The time jitter of 
the waveforms with relatively large amplitudes was approxi- 
mately 342ns(the o value with a Gaussian fit). 

Further analysis can be performed with time-walk correc- 
tion, which was not implemented in the FPGA in the pro- 
posed scheme. A cubic polynomial, shown in Figure 5 with a 
red solid line, was used to fit the 2-dimension histogram. The 
arrival time histograms with and without time-walk correc- 
tion are shown in Figure 6. Good correction can be obtained 
with small waveforms, and the time jitter is approximately 
344 ns(the o value with a Gaussian fit), which is almost iden- 
tical the value obtained with large waveforms. 
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Fig. 5. The relationship between waveform amplitude and time ob- 
tained with LED method. Red line shows the fit result with a cubic 
polynomial, used to correct the time jitter. 


D. Energy 


In general, two quantities are used to extract the energy 
loss from the detector: the amplitude and total charge of the 
signal. These two quantities were measured using ADCs and 
charge-to-digital converters (QDCs) in a conventional DAQ 
system. Using a digital approach, these measurements are 
replaced by appropriate algorithms in the FPGA. However, 
the signal extracted from the CSA was too wide, indicating 
that the digital QDC method was not suitable for the CsI(TI) 
array. 
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Fig. 6. (Color online) Time spectra obtained using LED technique. 
Blue dash line is the spectrum without time correction, red solid line 
is the spectrum with time correction. 


The most direct approach to extract the signal amplitude 
is to determine the maximum( positive signal) or minimum( 
negative signal) waveform. However, the measured ampli- 
tudes were significantly influenced by noise. To improve 
the energy resolution, digital filters, which perform the same 
function as shaping amplifiers in an analog measurement sys- 
tem, were used to shape the signal before amplitude extrac- 
tion. Theoretically, the best signal shape that maximizes the 
SNR is the infinite-width cusp [42]. However, a practical fil- 
ter of this type is the finite-width cusp filter that limits the 
amplitude measurement of a single signal to a specific time. 
The algorithm uses a different function instead of an infinite 
exponential function and performs truncation [43, 44], which 
means that the performance is reduced. Other commonly used 
shaping filters include a series of trapezoidal [45-51] and 
CR-RC”™ filters [52-54]. To identify a suitable shaping fil- 
ter, the y-rays produced by a ©°Co source were measured on 
a test bench operating in waveform mode, and the recorded 
data were processed with various digital filters. The wave- 
forms before and after applying the shaping filters are shown 
in Figure 7. All parameters in each filter were optimized us- 
ing repeated trials. The energy resolutions of the full-energy 
peaks for 1.33 MeV yrays were calculated as the criterion, 
and the results are shown in Table 2. The energy resolutions 
achieved by each filter were comparable with the performance 
of the finite cusp filter being slightly better. This is because 
the noise in the output test was extremely low, and the SNR 
values did not improve significantly. It should be noted that 
the noise level in the test is of the same order of magnitude 
as that of the ETF, indicating that all filters meet the require- 
ments in terms of performance. 

In [55], it is concluded that the family of trapezoidal 
filters offers good performance and simpler implementation 
among many other shaping filters, which influences this 
choice. Many algorithms exist for the implementation of 
trapezoidal filters, of which moving-window deconvolution 
(MWD) [49-51] is the simplest for FPGA implementation. 
Therefore, the MWD filter was selected as the shaping fil- 
ter for the CsI(TI) array, and the corresponding °°Co energy 
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spectrum is indicated by the red dashed line in Figure 12. 


E. PID algorithms 


As previously mentioned, performing a pulse shape analy- 
sis for the CsI(TI) array can reduce the background of charged 
particles in the energy spectrum of y-rays and improve the 
PID performance through a combination with AE-E method. 
The basis of pulse shape analysis is that the ratio of the fast 
and slow components of light generated by the CsI(T]) crys- 
tal depends on the type of incident particles, resulting in dif- 
ferent shapes of the output waveform. Multiple methods ex- 
ist for extracting these differences, among which the digital 
charge comparison [56-59] and rise-time comparison meth- 
ods [59, 60] are widely used and straightforward to imple- 
ment in the FPGA. Another method worth focusing on is 
called reconstructive particle identification (RPID) [51]. One 
reason is that this method was developed for CALIFA, whose 
construction is similar to the CsI(T1) array. The RPID method 
can be successfully migrated into the DAQ system with sig- 
nificant potential. Another reason is that RPID can directly 
extract the fast and slow components of the CsI(T1) crystal, 
thereby achieving a more improved performance. 

To compare these PID algorithms, a triple œ source 
(244Cm, 78°Pu, and 241Am) and a ©°Co source were used 
separately to irradiate the elemental detector of the CsI(TI) 
array under identical conditions. The CsI(T1) crystal enve- 
lope was pierced with a small hole to allow the a-particles to 
access the crystal. With respect to the digital charge compar- 
ison method, because the input signals are extracted from the 
CSA, which results in information on the incident particles 
being contained in the leading edges of the digitized wave- 
forms, both time windows for charge integration were set to 
include the leading edge of the waveform with the same start- 
ing point, as shown in Figure 8. The starting point is deter- 
mined using the dCFD method without interpolation. The in- 
tegration windows size was determined through repeated tri- 
als. With respect to the rise-time comparison method, the rise 
time is defined as the time frame between the points where 
the waveforms cross 12.5 and 87.5% of the full amplitude in 
the leading edges. Because the method performance is cor- 
related with the time accuracy, the linear interpolation tech- 
nique was used to determine the cross-points. Regarding the 
RPID method, some parameters are identical to those of the 
MWD filter previously described, with three additional pa- 
rameters: the fast and slow decay times of the CsI(TI) crystal 
and the time window for the second MWD procedure. The 
fast and slow decay times are in accordance with [51], be- 
cause they are almost constants. The length of the second 
time window was determined using repeated trials. 

The two-dimensional PID spectra for each method are 
shown in Figure 9. Superficially, œ and y-rays were well 
identified for all three methods. The spectrum shown in 
Figure 9(a) illustrates that the rise time and the amplitude 
of the waveform are weakly correlated and can be consid- 
ered independent. The red dashed lines in Figures. 9(b) 
and (c) indicate that the spectra approach zero when both 
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the horizontal and vertical axes are reduced. These rela- 
tionships indicate that the PID parameters listed in Table 3 
can be used to convert two-dimensional spectra into one- 
dimensional histograms, and the figure-of-merit (FoM), as 
defined in [61, 62], can be used to quantify the separation 
performance of the PID methods. These one-dimensional his- 
tograms are shown in Figure 10 and the relevant FoMs are il- 
lustrated. Although data in the low-energy regions, which can 
make the FoM values larger than those in the beam experi- 
ments, were not included, the FoM values shown in Figure 10 
remain good references for evaluating these three methods. 
Of the three methods, the rise time comparison achieved the 
worst score, which is consistent with the results in [28]. This 
is because the accuracy of the rise time deteriorated as the am- 
plitude of the input signal decreased. Surprisingly, the RPID 
method does not work as well as the digital charge compari- 
son method. One reason for this is the approximate treatments 
in the algorithm implementation process, such as calculating 
the exponential function. Moreover, it was identified that the 
decay time of the fast component in the CsI(T1) array depends 
on the type of incident particles (more precisely, on the aver- 
age energy density deposited in the crystal) [22, 63] and even 
on the total energy of the particles [20, 21], while the slow 
component decay time does not. If this is the case, RPID may 
not be an accurate method because of the underlying assump- 
tion of constant decay time. However, this remains sufficient 
for the separation of y-rays and light-charged particles. 

Overall, the aforementioned descriptions support the con- 
clusion that the digital charge comparison method is preferred 
owing to its good performance and ease of FPGA implemen- 
tation. Performance of the RPID method lags slightly; how- 
ever, the complexity of this algorithm rules this option out. 
The rise-time comparison method had the worst performance 
among the three methods; however, the result remains accept- 
able, making it an alternative option to the digital charge com- 
parison method. 


IV. FINAL SCHEME AND PERFORMANCES OF DPP 
ALGORITHMS IN THE FPGA 


Considering these aforementioned points, the final DPP al- 
gorithm scheme which includes a moving average filter, base- 
line restoration, leading-edge discrimination, moving win- 
dow deconvolution, and digital charge comparison, is formed. 
Thus, the diagram in Figure 1 can be improved as shown in 
Figure 11. All DPP algorithms are "FPGA friendly,” and fur- 
ther simplification can be achieved for the MAF and MWD 
by changing the algorithms to the recursive form. The algo- 
rithms were implemented in the DAQ system FPGA, com- 
pleting the algorithm mode. The sampling rate was reset to 
50 MS/s because the amount of data is significantly reduced. 
Subsequently, all procedures in Figure 11 are retested in this 
mode with the same radioactive sources used in Section IIE. 
The energy spectra are shown in Figures. 12 and 13. A 
small shift can be observed between the two energy spectra 
of the ©°Co source compared to the result obtained in the raw 
waveform mode. This is because of the rounding operation 
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TABLE 2. The performances of different shaping filters. 


Energy resolution 


Filters (FWHM, 1.33MeV) References 
Cusp 5.58% [43] 
Trapezoidal (single exponential function) 5.74% [47] 
Trapezoidal (bi-exponential function) 5.85% [47] 
Trapezoidal (MWD) 5.67% [51] 
CRpz-RC2 5.70% [53] 
CRpz-RC4 5.79% [53] 
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Fig. 7. (Color online) Shapes of raw waveform (a) and those modified by the following filters: (b) cusp, (c) trapezoidal with a single 
exponential function, (d) trapezoidal with a bi-exponential function, (e) MWD, (f) CRpz-RC2, and (g) CRpz-RC4. 


TABLE 3. The PID functions for three PID methods. 
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1.33 MeV full energy peak are almost unchanged. 


The performance of the PID controller with the on-line al- 
gorithms is shown in Figure 14. Compared with the same 
results obtained from the raw waveform mode shown in Fig- 
ure 10, almost the same positions where the peaks of the PID 
parameters are located can be found using the rise time and 
charge comparison methods. There is a slight improvement 
in the FoM parameter when comparing the charge compari- 
son methods. This is owing to an improvement in the DAQ 
system sampling frequency. For the rise time method, the 
FoM parameter shows little differences. The reason for this 
is the rounding operation for the final result in the FPGA. 
Therefore, it can be concluded that only a minimal perfor- 
mance difference exists between the on-line and off-line DPP 
algorithms. 
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Table 4 lists the key performance metrics for the proposed 
algorithm. Good energy resolutions and PID performances 
are achieved, indicating that the on-line algorithms in the 
FPGA are well formed, and the final DPP algorithm scheme 
can adequately meets the requirements. 
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Fig. 9. 2-dimensions PID spectra for three methods: (a) rise-time comparison, (b) digital charge comparison, and (c) RPID. The QLongGate 
and QShortGate are the two integral values in the long and short gates as shown in Figure 8. The N, and Nọș are the relative amplitude values 
of the slow and fast decay components of the CsI(T1) crystal.Red dashed lines in (b) and (c) indicate that the spectra approach zero when both 
horizontal and vertical axis values are reduced. The rise-time comparison spectrum shows that rise time and amplitude of the waveform are 
weakly correlated and can be considered as independent. This indicates that the list of PID parameters in Table 3 is well defined owing to the 


nearly constant values of the PID parameters. 
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Fig. 10. Histograms of PID parameters including FoMs for three methods: (a) rise-time comparison, (b) digital charge comparison, and (c) 
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Fig. 11. Final DPP algorithm scheme in FPGA for the CsI(T1I) array 
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TABLE 4. Major parameters obtained from Figures 12, 13 and 14 


by on-line algorithms. 


Parameters 

energy resolution(1.17 MeV gamma-rays) 
energy resolution(1.33 MeV gamma-rays) 
energy resolution(5.16 MeV a °3°Pu) 
energy resolution(5.48 MeV a 741 Am) 
energy resolution(5.80 MeV a *“*Cm) 
FoM(digital charge comparison) 
FoM(Rising time comparison) 


6.3% 
5.7% 
10.4% 
7.8% 
6.1% 
4.995 
2.329 
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Fig. 12. (Color online) Energy spectra of ®°Co radioactive source. 
The histogram shown by the blue solid line is obtained by the DAQ 
system algorithm mode and the on-line MWD algorithm, while the 
histogram shown by the red dashed line is obtained by the raw wave- 
form mode and the off-line MWD algorithm. 


V. SUMMARY 


In this study, a scheme for DPP algorithms was devel- 
oped for the CsI(Tl) array at ETF. A test bench with a and 
y sources was constructed to determine the algorithms in 
each step, resulting in the following final scheme: moving 
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Fig. 13. (Color online) Triple a sources energy spectra obtained by 
the on-line MWD algorithm. A triple Gauss function, shown by the 
red line, is used to fit the peaks. 


[1] 


[2] 


[3] 
[4] 


[5] 


[6 


=e 


[7] 


[8 


finer 


[9] 


[10] 


[11] 


[12] 


Isao Tanihata, Radioactive beam science, past, present and fu- 
ture. Nucl. Instrum. Methods Phys. Res. B 266, 4067 (2008). 
doi: 10.1016/j.nimb.2008.05.088 

Alexandra Gade, Nuclear spectroscopy with fast exotic 
beams. Phys. Scr. T152, 014004 (2013). doi: 10.1088/0031- 
8949/2013/T152/014004 

Y.Z. Sun, Study of Single Proton Knockout from ‘°C. Ph.D. 
Thesis, University of Chinese Academy of Sciences (2019). 
Y.Z. Sun, Z.Y. Sun, S.T. Wang et al., The charged frag- 
ment detector system of the External Target Facility. Nucl. 
Instrum. Methods Phys. Res. A 927, 390 (2019). doi: 
10.1016/j.nima.2019.02.067 

X.H. Zhang, S.W. Tang, P. Ma et al, A multiple sam- 
pling ionization chamber for the External Target Facility. 
Nucl. Instrum. Methods Phys. Res. A 795, 389 (2015). doi: 
10.1016/j.nima.2015.06.022 

Y. Sun, Z.Y. Sun, Y.H. Yu et al., Design and construction 
of a time-of-flight wall detector at External Target Facility of 
HIRFL-CSR. Nucl. Instrum. Methods Phys. Res. A 893, 68 
(2018). doi: 10.1016/j.nima.2018.03.030 

Y.Z. Sun, Z.Y. Sun, S.T. Wang et al., The drift cham- 
ber array at the External Target Facility in HIRFL-CSR. 
Nucl. Instrum. Methods Phys. Res. A 894, 72 (2018). doi: 
10.1016/j.nima.2018.03.044 

K. Yue, Z.Y. Sun, S.T. Wang et al., A CsI(Tl) gamma de- 
tection array at the external target hall of CSRm. Nucl. 
Instrum. Methods Phys. Res. B 317, 653 (2013). doi: 
10.1016/j.nimb.2013.07.037 

S. Takeuchi, T. Motobayashi, Y. Togano et al., DALI2: A 
Nal(Tl) detector array for measurements of y rays from fast 
nuclei. Nucl. Instrum. Methods Phys. Res. A 763, 596 (2014). 
doi: 10.1016/j.nima.2014.06.087 

Pieter Doornenbal, In-beam gamma-ray spectroscopy at the 
RIBF. Prog. Theor. Exp. Phys. 2012, 03C004 (2012). doi: 
10.1093/ptep/pts076 

H. Alvarez-Pol, J. Benlliure, E. Casarejos et al., Design 
studies and first crystal tests for the R°B calorimeter. Nucl. 
Instrum. Methods Phys. Res. B 266, 4616 (2008). doi: 
10.1016/j.nimb.2008.05.113 

H. Alvarez-Pol, N. Ashwood, T. Aumann et al., Performance 
analysis for the CALIFA Barrel calorimeter of the RB exper- 
iment. Nucl. Instrum. Methods Phys. Res. A 767, 453 (2014). 
doi: 10.1016/j.nima.2014.09.018 


4 


© 


1 
492 
493 
494 
495 
496 
497 
498 


average filter, baseline restoration, leading-edge discrimina- 
tion, moving window deconvolution, and digital charge com- 
parison. Subsequently, the DAQ system algorithm mode is 


completed using these DPP algorithms. 


It was identified 


that the performance does not change significantly between 
on-line and off-line algorithms. With the algorithm mode, 
good performances in the energy spectrum and PID, as listed 
in Table 4, are achieved, which indicates that the proposed 
49 DPP algorithm scheme meets the requirements to upgrade the 
soo CsI(Tl)array DAQ system at the ETF of the HIRFL-CSR. 


s44 [13] Lei Zhao, Long-Fei Kang, Jia-Wen Zhou et al, A 16- 


545 
546 
547 
548 
549 
550 
551 
552 
553 
554 
555 
556 
557 
558 
559 
560 
561 
562 
563 
564 
565 
566 
567 
568 


[14] 


[15] 


[16] 


[17] 


[18] 


Channel high-resolution time and charge measurement mod- 
ule for the external target experiment in the CSR of HIRFL. 
Nucl. Sci. Tech. 25, 010401 (2014). doi: 10.13538/j.1001- 
8042/nst.25.010401 

Xing-Wen Zhao, Yi Qian, Jie Kong et al., Readout elec- 
tronics for CSR-ETF silicon strip array detector system. 
Nucl. Sci. Tech. 25, 040402 (2014). doi: 10.13538/j.1001- 
8042/nst.25.040402 

Hai-Bo Yang, Xian-Qin Li, Yu-Hong Yu et al., Design and 
evaluation of prototype readout electronics for nuclide detec- 
tor in detector in Very Large Area Space Telescope. Nucl. Sci. 
Tech. 33, 65 (2022). doi: 10.1007/s41365-022-01047-5 
Yu-Ying Li, Chang-Yu Li, Kun Hu, Design and develop- 
ment of multi-channel front end electronics based on dual- 
polarity charge-to-digital converter for SiPM detector applica- 
tions. Nucl. Sci. Tech. 34, 18 (2023). doi: 10.1007/s41365-023- 
01168-5 

J. Kong, Y. Qian, H. Zhao et al., Development of the read- 
out electronics for the HIRFL-CSR array detectors. Journal 
of Instrumentation 14, P02012 (2019). doi: 10.1088/1748- 
0221/14/02/P02012 

Steven W. Smith, The Scientist and Engineer’s Guide to Digital 
Signal Processing, Second Edition. California Technical Pub- 
lishing, 1999. 


seo [19] Hongna Liu, Study of Neutron-Proton Correlations and 3N 


570 
571 
572 
573 
574 
575 
576 
577 
578 
579 
580 
581 
582 
583 
584 
585 
586 


[20] 


[21] 


[22] 


[23] 


Forces in !2C. Ph.D. Thesis, Peking University (2015). 
Agnieszka Syntfeld-Kazuch, Marek Moszyński, Lukasz 
Swiderski et al., Light Pulse Shape Dependence on y-ray En- 
ergy in CsI(Tl). IEEE Trans. Nucl. Sci. 55, 1246 (2008). doi: 
10.1109/TNS.2008.922805 

X. Lu, S. Gridin, R.T. Williams et al., Energy-Dependent Scin- 
tillation Pulse Shape and Proportionality of Decay Compo- 
nents for CsI:Tl: Modeling with Transport and Rate Equations. 
Phys. Rev. Appl. 7, 014007 (2017). doi: 10.1103/PhysRevAp- 
plied.7.014007 

R.S. Storey, W. Jack, A. Ward, The Fluorescent Decay of 
CsI(Tl) for Particles of Different Ionization Density. Proc. 
Phys. Soc. 72, 1 (1958). doi: 10.1088/0370-1328/72/1/302 

H. Grassmann, E. Lorenz, H.-G. Moser, Properties of CsI(T1) 
- Renaissance of an old scintillation material. Nucl. Instrum. 
Methods Phys. Res. A 228, 323 (1985). doi:10.1016/0168- 
9002(85)90276- 1 


Amplitude (a.u.) 


1 Ni al 1 
2000 2500 3000 3500 


Risetime (ns) 


f 1 
1000 1500 4000 


© FoM = 2.322 


Ol 


Y 


| ni 1 
1500 2000 2500 3000 3500 


PID parameter (rise time / ns) 


1000 4000 


5837 [24] P. Schotanus, R. Kamermans, P. Dorenbos, Scintillation char- 
588 acteristics of pure and Tl-doped CsI crystals. IEEE Trans Nucl. 
Sci. 37, 177 (1990). doi: 10.1109/23.106614 
S. Carboni, S. Barlini, L. Bardelli et al., Particle iden- 
tification using the AE-E technique and pulse shape dis- 
crimination with the silicon detectors of the FAZIA project. 
Nucl. Instrum. Methods Phys. Res. A 664, 251 (2012). 
doi: 10.1016/j.nima.2011.10.061 
J. Alarja, A. Dauchy, A. Giorni et al., Charged parti- 
cles identification with a CsI(T1) scintillator. Nucl. Instrum. 
Methods Phys. Res. A 242, 352 (1986). doi:10.1016/0168- 
9002(86)90232-9 
S. Aiello, A. Anzalone, G. Cardella et al., Light response and 
particle identification with large CsI(Tl) crystals coupled to 
photodiodes. Nucl. Instrum. Methods Phys. Res. A 369, 50 
(1996). doi:10.1016/0168-9002(95)00763-6 
[28] W. Skulski, M. Momayezi, Particle identification in CsI(T1) us- 
ing digital pulse shape analysis. Nucl. Instrum. Methods Phys. 
Res. A 458, 759 (2001). doi: 10.1016/S0168-9002(00)00938-4 
sos [29] Abraham Savitzky, Marcel J.E. Golay, Smoothing and Dif- 
607 ferentiation of Data by Simplified Least Squares Procedures. 
608 Anal. Chem. 36, 1627 (1964). doi:10.1021/ac60214a047 
sos [30] P. Marchand, L. Marmet, Binomial smoothing filter: A way 
610 to avoid some pitfalls of least-squares polynomial smoothing. 
611 Rev. Sci. Instrum. 54, 1034 (1983). doi:10.1063/1.1137498 
612 [31] Paul H.C. Eilers, A Perfect Smoother. Anal. Chem. 75, 3631 
613 (2003). doi:10.1021/ac034173t7 
614 [32] É. Barat, T. Dautremer, T. Montagu, et al, A bi- 
615 modal Kalman smoother for nuclear spectrometry. Nucl. 
616 Instrum. Methods Phys. Res. A 567, 350 (2006). doi: 
61 10.1016/j.nima.2006.05.243 
618 [33] David G. Abrecht, Jon M. Schwantes, Ravi K. Kukkadapu, 
619 et al., Real-time noise reduction for Méssbauer spectroscopy 
620 through online implementation of a modified Kalman filter. 
621 Nucl. Instrum. Methods Phys. Res. A 773, 66 (2015). doi: 


[25] 


[26] 


[27] 


605 


rá 


622 
623 
624 
625 
626 
627 
628 
629 
630 
631 
632 
633 
634 
635 
636 
637 
638 
639 
640 
641 
642 
643 
644 
645 
646 
647 
648 
649 
650 
651 
652 
653 
654 
655 


10 


10° 


(b) 


1 1 
10000 40000 50000 60000 


[34] 


[35] 


[36 


ol 


[37] 


[38 


peer 


[39] 


[40] 


[41] 


(d) 


0.08 


0.1 0.12 0.14 
PID parameter(Q_/Q, ) 


Fig. 14. On-line PID algorithm performance of digital charge comparison and rise-time comparison method in FPGA. 


10.1016/j.nima.2014.10.053 

A. Geraci, I. Rech, E. Gatti, et al., Shared baseline restora- 
tion at minimum noise for high resolution spectroscopy. 
Nucl. Instrum. Methods Phys. Res. A 482, 441 (2002). doi: 
10.1016/S0168-9002(01)01509-1 

E.M. Khilkevitch, A.E. Shevelev, I.N. Chugunov, et al., 
Advanced algorithms for signal processing scintillation 
gamma ray detectors at high counting rates. Nucl. In- 
strum. Methods Phys. Res. A 977, 164309 (2020). doi: 
10.1016/j.nima.2020. 164309 

Jia-wen Zhou, Shu-bin Liu, Long-fei Kang, et al., Research 
on the Trigger System of Readout Electronics in HIRFL-CSR. 
Nuclear Electronics & Detection Technology 33, 152 (2013). 
Mark A. Nelson, Brian D. Rooney, Derek R. Dinwiddie, et 
al., Analysis of digital timing methods with BaF» scintillators. 
Nucl. Instrum. Methods Phys. Res. A 505, 324 (2003). doi: 
10.1016/S0168-9002(03)01078-7 

A. Fallu-Labruyere, H. Tan, W. Henning, et al., Time res- 
olution studies using digital constant fraction discrimination. 
Nucl. Instrum. Methods Phys. Res. A 579, 247 (2007). doi: 
10.1016/j.nima.2007.04.048 

L. Bardelli, G. Poggi, M. Bini, et al., Time measurements by 
means of digital sampling techniques: a study case of 100 ps 
FWHM time resolution with a 100 MSample/s, 12 bit digitizer. 
Nucl. Instrum. Methods Phys. Res. A 521, 480 (2004). doi: 
10.1016/j.nima.2003.10.106 

K. Wang, S. Samaranayake, A. Estrade. Investigation of a digi- 
tizer for the plastic scintillation detectors of time-of-flight mass 
measurements. Nucl. Instrum. Methods Phys. Res. A 1027, 
166050 (2022). doi: 10.1016/j.nima.2021.166050 

Shefali Saxena, and Ayman I. Hawari, Investigation of FPGA- 
Based Real-Time Adaptive Digital Pulse Shaping for High- 
Count-Rate Applications. IEEE Trans. Nucl. Sci. 64, 1733 
(2017). doi: 10.1109/TNS.2017.2692219 


~ 682 


656 [42] V. Radeka, Optimum Signal-Processing for Pulse-Amplitude 
Spectrometry in the Presence of High-Rate Effects and 
Noise. IEEE Trans. Nucl. Sci. 15, 455 (1968). doi: 
10.1109/TNS.1968.4324970 
Yunchen Qian, Huaiqiang Zhang, Fenhua Lu, et al., Parameter 
optimization and pile-up identification of cusp shaping for nu- 
clear pulse signal. Nuclear Techniques 44(11), 110402 (2021). 
doi: 10.11889/j.0253-3219.2021.hjs.44.110402 
M. Bogovac and C.Csato, Implementation of a truncated cusp 
filter for real-time digital pulse processing in nuclear spectrom- 
etry. Nucl. Instrum. Methods Phys. Res. A 694, 101 (2012). 
doi: 10.1016/j.nima.2012.07.042 
V. Radeka, Trapezoidal filter of signals from large germanium 
detectors at high rates. Nucl. Instrum. Methods 99, 525 (1972). 
doi: https://doi.org/10.1016/0029-554.X(72)90666-0 
Valentin T. Jordanov and Glenn F. Knoll, Digital synthe- 
sis of pulse shapes in real time for high resolution radiation 
spectroscopy. Nucl. Instrum. Methods Phys. Res. A 345, 337 
(1994). doi: 10.1016/0168-9002(94)91011-1 
[47] Cosimo Imperiale and Alessio Imperiale, On nuclear spectrom- 
6 etry pulses digital shaping and processing. Measurement 30, 49 
(2001). doi: 10.1016/S0263-2241(00)00057-9 
3 [48] Huai-Qiang Zhang, Lian-Quan Ge, Bin Tang et al., Optimal 
679 choice of trapezoidal shaping parameters in digital nuclear 
680 spectrometer system. Nucl. Sci. Tech. 24, 060407 (2013). doi: 
10.13538/j.1001-8042/nst.2013.06.011 
[49] Andrey Georgiev and Werner Gast, Digital Pulse Processing in 
683 High Resolution, High Throughput Gamma-Ray Spectroscopy. 
684 IEEE Trans. Nucl. Sci. 40, 770 (1993). doi: 10.1109/23.256659 
685 [50] J. Stein, F. Scheuer, W. Gast, et al., X-ray detectors with digi- 
686 tized preamplifiers. Nucl. Instrum. Methods Phys. Res. B 113, 
687 141 (1996). doi:10.1016/0168-583X(95)01417-9 
ess [51] M. Bendel, R. Gernhäuser, W.F. Henning, et al., RPID - A new 
689 digital particle identification algorithm for CsI(T]) scintillators. 
690 Eur. Phys. J. A 49, 69 (2013). doi: 10.1140/epja/i2013-13069-8 
691 [52] M. Nakhostin, Recursive Algorithms for Real-Time Digital 
692 CR-(RC)” Pulse Shaping. IEEE Trans. Nucl. Sci. 58, 2378 
693 (2011). doi: 10.1109/TNS.2011.2164556 
« 694 [53] Yinyu Liu, Jinglong Zhang, Lifang Liu, et al., Implementa- 
695 tion of real-time digital CR-RC™ shaping filter on FPGA for 
696 gamma-ray spectroscopy. Nucl. Instrum. Methods Phys. Res. 


657 
658 
659 
e60 [43] 
661 
662 
663 
664 [44] 
665 
666 
667 
668 
669 
670 


[45] 


671 [46] 
672 
673 
674 
675 


| E 


N 


677 
6 


N 


681 


697 
698 
699 
700 
701 
702 
703 
704 
705 
706 
707 
708 
709 
710 
711 
712 
713 
714 
715 
716 
717 
718 
719 
720 
721 
722 
723 
724 
725 
726 
727 
728 
729 
730 
731 
732 
733 
734 
735 
736 
737 


[54] 


[55] 


[56] 


[57] 


[58] 


[59] 


[60] 


[61] 


11 


A 906, 1 (2018). doi: 10.1016/j.nima.2018.05.020 
Huai-Qiang Zhang, Zhuo-Dai Li, Bin Tang, et al., Optimal pa- 
rameter choice of CR-RC” digital filter in nuclear pulse pro- 
cessing. Nucl. Sci. Tech. 30, 108 (2019). doi: 10.1007/s41365- 
019-0638-7 

J. Kamleitner, S. Coda, S. Gnesin, et al., Comparative anal- 
ysis of digital pulse processing methods at high count rates. 
Nucl. Instrum. Methods Phys. Res. A 736, 88 (2014). doi: 
10.1016/j.nima.2013.10.023 

M. Moszynski, D. Wolski, T. Ludziejewski, et al., Parti- 
cle identification by digital charge comparison method ap- 
plied to CsI(TI) crystal coupled to photodiode. Nucl. Instrum. 
Methods Phys. Res. A 336, 587 (1993). doi: 10.1016/0168- 
9002(93)91267-Q 

Y. Kaschuck and B. Esposito, Neutron/y-ray digital pulse 
shape discrimination with organic scintillators. Nucl. In- 
strum. Methods Phys. Res. A 551, 420 (2005). doi: 
10.1016/j.nima.2005.05.071 

Jamil-Qureshi Faisal, Jian-Ling Lou, Zhi-Huan Li et al., 
A pulse shape discrimination of CsI(TI) crystal with °He 
beam. Nucl. Sci. Tech. 21, 35 (2010). doi: 10.13538/j.1001- 
8042/nst.21.35-38 

Chang-Lin Lan, Xi-Chao Ruan, Gang Liu et al., Particle 
identification using CsI(T1) crystal with three different meth- 
ods. Nucl. Sci. Tech. 19, 354 (2008). doi: 10.1016/S1001- 
8042(09)60018-X 

S.D. Jastaniah and P.J. Sellin, Digital techniques for n/y pulse 
shape discrimination and capture-gated neutron spectroscopy 
using liquid scintillators. Nucl. Instrum. Methods Phys. Res. A 
517, 202 (2004). doi: 10.1016/j.nima.2003.08.178 

R.A. Winyard, J.E. Lutkin and G.W. McBeth, Pulse shape 
discrimination in inorganic and organic scintillators,I. Nucl. 
Instrum. Methods 95, 141 (1971). doi: 10.1016/0029- 
554X(71)90054-1 

R.A. Winyard and G.W. McBeth, Pulse shape discrimination in 
inorganic and organic scintillators,II. Nucl. Instrum. Methods 
98, 525 (1972). doi: 10.1016/0029-554X(72)90238-8 

J.A. Biggerstaff, R.L. Becker and M.T. Mcellistrem, Charged 
particle discrimination in a CsI(T1) detector. Nucl. Instrum. 
Methods 10, 327 (1961). doi: 10.1016/S0029-554X(61)80127- 
4 


